On a Local Concept of Wave Velocities 



I. V. Drozdov * 

and 

A. A. Stahlhofen t 

University of Koblenz, Institute for Natural Sciences 
Department of Physics 
Universitatsstr.l, D-56080 Koblenz, Germany 

February 7, 2008 



Abstract 

The classical far field concept of wave velocities has its merits while 
exibiting intrinsic difficulties. 

A general local approach for the definition of velocities and especially 
phase velocities for waves avoiding these difficulties is proposed. It includes 
the classical definitions as particular cases and can be applied to waves of 
an arbitrary structure, and to arbitrary propagation media as well. Appli- 
cations of the formalism are elucidated and some basic properties of the 
local concept defined here are discussed. 



Some remarks on wave velocities 



Waves are conventionally described via propagating harmonic functions contain- 
ing a periodic factor of the form e lLp . The argument of these periodic functions , 
interpreted as the "phase", is usually identified with (t — x/c) via the free wave 
equation, where c is the "phase velocity" []]. 

This canonical approach originates from classical wave optics and analyzes 
ordinary light waves propagating in the "far field", where the wave is far away 
from the source. Such a light wave is a solution of a source free wave equation 
described by the periodic function mentioned above. The Ansatz earned his own 
merits in the context of several wave phenomena and in particular classical wave 
optics [2]. 
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A second fundamental concept of characterizing wave propagation is that of 
the "group velocity". This concept is - mathematically speaking - a comprehen- 
sive definition for a specified (linear) superposition of solutions of the free wave 
equation with the same periodicity properties usually expressed by the frequency 
distribution of the constituent periodic waves. This definition is also an inher- 
ent far field concept considering "source-free" waves propagating in a strongly 
homogeneous and isotropic medium. This medium is characterised only by its 
"dispersion" (supposed to satisfy additionally the Kramers- Kronig relations), i.e. 
by a dependence of the complex index of refraction n(oS) on the periodicity pa- 
rameter of the wave - the frequency u. 

This basis for definitions of "phase" and "group" velocities turns out to be 
appropriate only for the special class of wave propagation mentioned above. The 
notions of a signal and its velocity developed in this context did accomodate 
many experimental data (see e.g. |6J). It is however by no means obvious, to what 
extent these far field concepts can be applied to near field problems. This leads 
to the question, if this classical approach is flexible enough to cover more general 
wave phenomena apart from the "special cases" mentioned above. A possible 
answer to this question is the focus of the present work. 

Before proceeding we list some actual problems accompanying the classical defi- 
nitions and suggesting a resolution by a local definition of a propagation velocity 
[T2] as proposed here. 

Theories developed for general not necessary periodic electromagnetic pulses 
[3[ [10], did draw basically on the adopted "canonical" scenario. The analysis 
for a wave of an arbitrary form, for instance, has been based in general (by 
analogy to the classical case) on a representation via periodic functions by means 
of Fourier analysis. A well defined "phase velocity" is assigned to each Fourier 
component equipped with its frequency; a corresponding "group velocity" has 
subsequently been defined for the complete Fourier superposition (i.e. for the 
"wave packet"). As a consequence, the classical wave theory (in all its facets) 
is based on apparently "canonical" definitions relying on the phase of periodical 
functions and of groups of such functions. Such a strategy treates this as the 
inherent kernel of wave phenomena. 

A close inspection, however, reveals mathematical deficiencies and several 
shortcomings with respect to physics thereby rendering the classical "canonical" 
approach to some extent non-natural as elucidated now. 

First of all, a propagating wave is per definition a local space-time distribution 
which is a solution of a local differential equation. Hence the definition of velocity 
as a space-time relation should respect this frame and can be expected to be local 
as well. 

The canonical definitions of a velocity ignore this fact and involve a frequency 
and a wave vector, which are nonlocal parameters in following the space- and 
time-periodicity, assumed separately in each case. Moreover, these parameters 
have nothing to do with the propagation process itself, since they do not enter 
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in the wave equation. Basically, the wave equation admits a propagation of ar- 
bitrary shapes without presumption of periodicity. In other words, the canonical 
approach bears anyway on a representation of an arbitrary solution by a set of 
periodic harmonic functions, each of that does not obey, generally speaking, the 
original wave equation (i.e. is not a solution). An attempt to apply it on waves re- 
sulting from non-linear equations such as solitons for instance, demonstrates this 
problem in a clear way since in this case even a linear superposition of solutions 
is not anymore a solution! 

Second, the subject of transmission and velocity of a signal is based on suf- 
ficiently local procedures of measurement [I]. Roughly expressed, the intervals 
between several space-time points are measured, where the attribute under study 
is detected. Definitions based on periodicity parameters can possess generically 
neither time- nor space-locality. Moreover, any Fourier transformation is basically 
a global object, and all manipulations concerned are in general mathematically 
exact only with integration over the whole space and whole time, as well as over 
an infinite frequency band for backward transformations. 

The "canonical" definitions of the "phase" and "group" velocity mentioned 
above pQ, are based on the special case of a propagation of periodic harmonic 
waves in homogeneous media. Thus there is already an essential contradiction 
between the local character of wave propagation and the inherent globality con- 
tained in the definitions of wave velocities. 

Any attempt to construct a local measurable object using Fourier sums or 
integrals contain an essential contradiction as outlined above and provides indeed 
no real locality. This is the source of several problems arising when replacing 
originally local features by " microglobal" ones [31 171 [TU]. 

The validity of this approach has to be checked from the mathematical point 
of view in the sense of functional analysis as well as from physical consistence 
and in fact it turns out to be correct only in special cases as mentioned above. 

Even if the approach is supported by mathematical consistency (like the as- 
sumption of an infinite frequency band pQ), it does not lend itself easily to a 
transparent physical interpretation an the calculations are cumbersome. 

These drawbacks, are for instance, the origin of all problems involved in the 
theory of signal transmission in media and an interpretation of the results. 

The classical papers since [I] and later improvements thereof [10J [H], still 
contain an essential mismatch between local and global wave features, based 
on several misleading definitions (like signal velocity, energy velocity and group 
velocity E]). This ansatz is bound to lead in applications to controversial 
results. For example, it seems not more to be surprising, that the subject of 
"group velocity" failes to describe a propagaton of ultra-short pulses [8]. 

These remarks shall suffice to motivate some alternative concept of wave prop- 
agation and wave velocities for the following reasons: 
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1. The canonical "harmonic" approach represents local wave features in terms 
of generically non-local attributes [121 [13]. Hence: 

2. The mathematical equipment is not exact for finite physical values. 

3. It is therefore difficult to impossible to apply it to for near-field effects, if 
the size of space-time regions are of an order of magnitude smaller than the 
wave periodicity parameter; 

4. It is poorly suited for inhomogeneous and anisotropic media, as well as for 
ultra-short pulses of arbitrary form and for nonlinear waves (7[ [15] . 

5. As a consequence, the applicabilty of this approach for several fields of 
modern quantum optics, nanooptics and photonics is barely justified, since 
these topics deal with the parameter areas outside its range of validity [9] . 

Moreover there is a clear reason to avoid this "harmonic" approach from the 
pure technical point of view: on the one hand it turns out to be an essential re- 
striction of generality and universality of the theory; on the other hand it brings 
to play ( especially the representation via Fourier series and integrals )a number 
of problematical artifacts such as an apparent violation of causality, infinite fre- 
quency bands for signals and so forth. These problems enforce further theoretical 
constructs like the analysis of dispersion relations [H] to resolve these artifical 
problems. 

A first step towards a more general concept should be based on an independent 
alternative approach without doubting the canonical "harmonic" criteria, leading 
to the same (verified) results in known special cases. 

The aim of the present paper is to establish a transparent criterion for evalu- 
ation of the propagation velocity for a quite arbitrary signal, that does not need 
an explicit representation in terms of harmonic or exponential functions at all, 
and with minimal loss of generality in other respects. First of all, we have to 
recall, what is being measured in experiments and what is meant when speaking 
about a "wave velocity", thereby providing the spectrum of wave velocities. This 
discussion is presented in Sec. 2. Sec. 3 is devoted to a test of the definition using 
some simple well known examples. Some remarkable properties of the behaviour 
of the wave velocities defined here under relativistic transformation are outlined 
in Sec. 4. The discussion is concluded in Sec. 5 by a practical interpretation and 
conversion from a local to a global evaluation of signal velocities. 

2 N-th order phase velocity (PV) 

The following discussion is restricted to a 1 + 1 dimensional space-time for sim- 
plicity. To evaluate the speed of a moving matter point one has to check the 
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change of the space coordinate Ax during the time interval At. How can this 
ansatz be applied to waves ? 

An arbitrary wave is a function defined on the 2-dimensional manifold (x, t) 
and one has no a priori defined fixed points to follow up as in the former case. 
Let us therefore analyze a continuous smooth function ip(x, t) of arbitrary shape. 

In a first step towards a local definition of velocity we consider a vicinity 
U(x) of a certain point x at the certain time t. Further we assume the local 
information about the function ip(x,t) to be measurable, i.e. we should be able 
(at least in principle) to evaluate the values of the function ip(x,t) itself and all 
its derivatives in the point x as well, and in the vicinity U{x). 

A description of propagation is based on monitoring the fate of a certain 
attribute ("labelled point") of this shape, i.e. one finds the same attribute at the 
next moment t\ in the next point X\. 
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Figure 1: The propagation of a signal ^(x, t) without a change of shape with the value 
ipo of the amplitude as a traced attribute. 

Let this attribute be labelled by some one-point fixed value tpo of the func- 
tion ip(x,t). Suppose, one follows up this local "attribute" and manages a local 
observation of the condition: 

^(x,t)-^ = 0, (2.1) 

which fixes the space-time points {x,t} where this condition holds (FigHJ). 
Thus we have a function x(t) given in an implicit form (12. ip . 
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The first order implicit derivation provides the velocity of this one-point local 
attribute: 

j (hp 

«<o)(M):=f = -§, (2-2) 

dx 

called from here on the "zero-order" or "one-point" phase velocity (0-PV). It de- 
scribes in the simplest case the speed of translation of some arbitrary pulse, that 
can be treated as the signal propagation velocity, provided that the measured 
value ipo of the amplitude is the signal considered. 

To proceed with the local description of the shape propagation we now con- 
sider as the measurable attribute of ip(x,t) at some time to not the single value 
if)o = ip(x ; to) at the point x , but a set of values ip in a local neighborhood 
(like U(x)) close to this point. Then another local attribute can be constructed 
to trace their propagation. If, for instance, one looks at a certain value of the 
first derivative dip(x;to)/dx := xo °f the shape ■?/>, (like di[)(x;to)/dx = 0, i.e. 
at a maximum or minimum point), one should consider the propagation of the 
condition: 

d^(x, t) 



Xo; (2.3) 



dx 

this provides the propagation velocity via 

dx 2 

called in view of (12.21) the "first order" or "two-point" phase velocity (1-PV) 
respectively. 

When tracking the propagation of a maximum (minimum), the conditions 
d^fc^ = and ^rt < have to hold simultaneously. 

ox ax z J 

This Ansatz is easily iterated to phase velocities of order N interpreted as the 
propagation velocities of higher order local shape attributes via 

d N + 1 j> 

leading to N-th order ("N+l-point") PV, describing the propagation of higher 
order local shape attributes. 

The phase velocities so defined are obviously local features depending on 
space-time coordinates. It has to be noted, that any given problem at hand 
might require a particular choice of a PV allowed for by the definitions given 
above. The following examples elucidate these requirements and demonstrate 
the flexibility of this definitions. 

Before proceeding to it should be recalled, that the PV-spectrum has been 
obtained wanting to describe the propagation of an arbitrary pulse in terms of 
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local attributes in a medium whose properties are dependent on several variables, 
especially time and space coordinates being the most prominent and natural 
examples thereof. Any initial shape ip(xo, to) thus should be deformed during the 
propagation (or evolution, as typicaly encountered for dispersive media). The 
ordinary phase velocity vq therefore is not a relevant criterion to characterize the 
shape propagation and one has to choose an appropriate PV v^y 

For example, let the pulse ip(x,t ) be subject to damping (FigfJ]). As a con- 
cequence, the zeroth order PV measured in the point X\ gives a magnitude much 
smaller as the same magnitude measured in the point x^- 

For an amplified signal (FigEJ like a signal propagating in a laser excited 
medium), by comparison, the zero order PV from the point X\ provides alto- 
gether even a backward propagation. In both cases an appropriate approach 
would be to apply the first order PV which describes the propagation of the 
maximum up from the point Xq properly. 




Figure 2: A damping deformation of a signal ip(x,t). A relevant propagation attribute 
is the peak location (maximum) 

For a propagation of a kink front that experiences a deformation, one can 
check the translation of the second derivative of the shape, keeping track of the 
turning-point of the kink shape (FigHj). In this case the second order ("three- 
point") PV turns out to be the relevant velocity of propagation. 

Finally, it should be noted, that the definiton of phase velocities of order zero 
and one, the and Vin respectively, admits a straightforward generalization 
to two-, three-, and higher-dimensional propagation, while the phase velocities of 
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Figure 3: The propagation of an amplified signal ip(x,t). The location of maximum 
describes the wave propagation 



second ( v nn) and higher orders inherently contain a certain element of ambiguity 
in definition, since a possibility to choose a second-order traced attribute is not 
unique [16J. 



3 Examples 



Let us consider the conventional 1+1- dimensional wave equation 

a 2 dt 2 dx 2 



ip(x,t) = 



possessing translational solutions of the form 



ip(x,t) = ip(t± 



(3.1) 



(3.2) 



It is easy to check that in this case the PV's of all orders defined above are 
identical and read 



v (N) (x,t) = a, AT = 0,1,2,... (3.3) 

which is nothing else but the classical phase velocity ±a, thereby satisfying 
the a priori definition of "phase velocity" itself as a medium constant in (13.1 1) . 

Let a propagating shape \l/ now be subject to a temporal damping similar to 
(Figj2]) with 



V(x,t) = i\){t - -)e~ M = M)e 
a 



S 




Figure 4: The propagation of a growing kink ("tsunami model"). The turning point 
is chosen to trace the propagation 



The ordinary 0-PV velocity reads 

w (0) =a(l- X^A , (3.5) 

where the prime denotes the derivative of ip with respect to its argument 
<p = t — x/a. 

This result provides a velocity with bad physical features: the velocity grows 
for a descending shape, for an ascending shape it decreases, can even be neg- 
ative, and it diverges exactly for the peak point (under the condition that the 
(measured) amplitude \I> as well as parameter a, X have positive values). 

A relevant physical velocity in this case is for instance the 1-PV 

v (I) =a(l- X^A (3.6) 

providing for a peak being traced exactly the canonical phase velocity. 
For this shape further PV's of higher orders are given by (|2.5p : 

v ( N) =a(l- X-^j = a (1 - A(log>W)- 1 ) (3.7) 

where ip( n > denotes the n-th derivative of ip with respect to its argument as 
mentioned above. For an amplified signal as in Figl3]we can e.g. change the sign 
of A. Especially, for the case of a kink (Fig. |4|) 

x 

ip(t ) = ip(4>) = arctan^, (3.8) 

a 

the spectrum of phase velocities reads by comparison : 



9 



f( ) = a (l + A(l + 2 ) arctan^)) , 

v (I) =a(l- >^f) , (3-9) 

The ordinary 0-PV has an oscillating sign at A and is multiple defined because 
of the arctan periodic function . Therefore it cannot be interpreted as a "well- 
defned" physical velocity. If the point being traced should be "labelled" by a 
derivative attribute, it turns out to be a physically inconvenient choice since the 
shape possesses no real peaks, their propagation could be traced. Moreover the 
velocity vm diverges at (ft = 0. 

The possible "labelled" attribute is also the turning-point traced by the 2-PV. 
The velocity vyi) also possesses two singularities at <p — ±l/\/3 which do not 
coincide with the labelled point <p = 0, so it can be successfully followed up at 
the measurement. 



Historical remark 



Canonically the definition of velocity should have proceeded with an originally 
artificial extension of a translational solution of ( 13. II) : 



tfj{x,t)=^(t--)=tp[-(ujt--x)], = <f)(ujt-kx), k = - (3.10) 
a to a a, 

and the parameter uj is further understood in restricted sense as a "frequency" 
of a "necessary" perodic harmonic function 0, usually e ±lx as mentioned above 
in the introduction. 

It is not surprising that the 0-order PV provides in this case the value ui/k, 
canonically interpreted as a "phase velocity of periodic wave" 

In the case of any adopted interrelations between k and uj that are not en- 
countered in the wave equation, in particular any so called "dispersion relation" 
between frequency and wave-number, the PV v (o) is in fact a proportionality 
factor 

aluj = V(tydk, (3-11) 
which is identical with the classical ("canonical") definiton of a "group veloc- 

ity" m. 

At this point it should be recalled, that the original idea of the "group veloc- 
ity" U, as summarized e.g. in [18J, 
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was of the similar form as the recent definition (12.21) of 0-PV; the equation (13.121) 
is however intrinsically controversal since being constructed of local derivatives of 
the non-local parameter A (wavelength). For a variable wavelength much smaller 
than a vicinity of the point, this "group velocity" provides therefore a natural 
approximation to the zero-order phase velocity t>(o). This was for a long time a 
reason to take it for a grainted relevant physical concept. 

In the present approach this feature appears as a physical phase velocity 
following in a straightforward way from the interpretation of phase propagation 
and does not require an interpretation of u and A; as a "frequency" and "wave 
number", as well as a constancy of some "group" respectively. 

Note, that in the case of kink there are no suitable definitions of a "wave- 
group" and of a concerned "group velocity" for this solution since a Fourier 
decomposition does not work on a non-compact support. 



4 Lorentz covariance 



The ordinary zero order PV (12.21 ) possesses a important property, i.e. a local 
covariance in the sense of special relativity, as shown below. This is not the case 
for PV's of higher orders. 

Actually the ordinary PV ivo) measured in some stationary system X takes 
in some other system X', moving with a constant speed V, via the Lorentz 
t r ansf or mat ions 

*' = 4^L, * = ( 4.i) 



the form 



c 2 

This means that the zero order PV respects the relativistic velocity addition. 
Especially, a subluminal zero order PV remains also subluminal in any other 
moving system X' . 

This result should not be a surprise, since the definition of the V(p) is implied 
by the condition: 

ip(x,t) = const, (4.3) 

which remains to be of the same form under arbitrary transormations x = 
x(x', t');t — t(x', t'), implying 

dip(x, t) = ^-dt + ^dx = (4.4) 
at ox 

for the first order differential form (or simply first differential), which possesses 
a form-invariance property under transformations. 
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Since 

V(o)(x,t) = — (4.5) 

per definition, the equation (14.41) inplies the definition (12.21) of 0-PV, so it should 
behave under space-time transformations as a usual velocity of a matter point. 

The first order PV (12 .4p . by comparison, evaluated in some system X trans- 
formes in the moving system X' to: 

1 + afst _ V \j?W + J?) } ^ 
J {i) ~ (v 2 a 2 , a 5 o v a 2 \ „,, ■ ^ ' 



If the pulse ip(x,t) obeys the free wave equation (I3.ip . the transformation 
becomes 



(l + £) v m + 2V 



" " + + <4 ' 7) 

i.e. a relation that should be called the "first order velocity addition". It dif- 
fers obviously from the corresponding transformation of Vm^, since the definition 
of the velocity Vm results from the condition 

<P**>)-0 (4,) 

whose form is explicitely non-invariant under space-time transformations. 

It can be shown that a subluminality of the first order PV is nevertheless still 
preserved by this transformation as well. Note, that for a signal which does not 
obey the free equation (13.11) . this restriction is in general not guaranteed anymore. 



5 A global velocity of signal transmission and 
"dynamic separation" 

The discussion of local velocities was aimed towards an evaluation of global fea- 
tures of signal propagation, namely the propagation through a finite spatial inter- 
val during a finite temporal interval. In other words, a global velocity, practically 
measured, means roughly the length of the distance Ax traveled by a traced at- 
tribute divided by the time interval At. 

The local PV's analyzed above can be interpreted as first order differential 
equations of the form 

V(N)(x,t) = —, (5.1) 
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Figure 5: Phase velocities as a family of isoclines V(m(x,t) and averaged global veloc- 
ities between two events (measurements) 



that can be illustrated graphically as a field of isoclines (FigfS]). 

Here the local PV is the tangent function of the tangent vector of the isocline, 
and the averaged (total) velocity between (t ,x ) and (ti,Xi) is represented by 
the tangent of the hypotenuse of the triangle {(t , x ), (t 1; Xi), (t , Xi)}, see Figj5j 

This procedure is elucidated best by some clear and well known examples: 
Consider the propagation of a translation mode of the form 

ip(x,t) = ip(£t - k(x)x) (5.2) 

describing, for instance, an electromagnetic wave propagation in an inho- 
mogeneous dielectric medium. The k(x) depends now on the space coordinate, 
resulting from the optical inhomogeneity: 

k(x) = (5.3) 

where n(x) is the spatially dependent index of refraction. Then the local zero 
order PV is provided by Eq. (12.21) and results in the first order equation 

>*><*.«> ^£ = -rV- = ^ (5 - 4) 

v ' at n'x + n ox 

An explicit integration of the (15.41) leads to 

c(ii — t ) = xin(xi) — xon(xo) (5.5) 
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for the interval of two events, where the signal attribute choosen for tracing 
has been checked at the time t in the point Xq and afterwards in x\ at t\. 
Assumed, the (to = 0,a;o = 0) is not a critical point (i.e. of the knot type) of 
the equation ( 15. ip . Then t\ = At, x\ = Ax and the global transition velocity 
between two points xq, x\ is 

v m {Ax) = (5.6) 

For the 1-PV velocity of eq. (12.4j) the same procedure leads to 

c c 

- log ((xn(x)) ) — (xn(x)) . (5.7) 



v (I) (x,t) £ 

This, equation straightforward to integrate, provides for the averaged velocity 
of two-point signal attribute along the distance Aa;: 

— = n(Ax) — \og(n'(Ax).Ax + n(Ax)) (5.8) 

v{i) £Ax 

It is worth to notice, that for the first order PV in media with a variable refraction 
index n(x) an essential dependence on the parameter £ enters. 

The meaning of £ can be derived from the given form of solution ip. It can be 
e.g. interpreted as a frequency factor for a periodic mode or a damping factor for 
evanescent one etc. This parameter simply means an enumeration of solutions of 
a solution family (space of solutions). The merit gained above is the separation of 
these solutions on the parameter £ through the different first (and higher)-order 
PV through the inhomogeneihty of media. 

In the case of a periodic wave for instance, it can be interpreted as a "dis- 
persion" although an explicit frequency-dependence of n(x) was not assumed. 
Moreover, the parameter £ has not necessarily to be interpreted as a "frequency" 
of some time periodic oscillation, but rather as a time component of the time- 
spatial wave vector {u;,k} for a translation mode (13.21) in a general form, thus 
the considered phenomenon has another origin and is much more general as a 
conventional non- localized frequency dispersion n(uS) 

Thus we established for the first (and higher) order PV in optical inhomoge- 
neous media the essential dependence on the t-component u of the wave transla- 
tion vector, especially frequency, even for local non dispersive media, which could 
be called " dynamic separation" . The global averaged dynamic separation of eq. 
(15. 8p survives as a corollary of the local separation of (15. 7p . This phenomena does 
not occur for the ordinary zero order PV. 



6 Concluding remarks 

The inherent inconsistency of the classical subjects of "phase velocity" and "group 
velocity", as well as "signal velocity" based therein has been discussed. The 
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inapplicability of these concepts for actual studies in photonics, near-field and 
nano-optics has been shown to result from the essential non-locality of these 
terms. 

An alternative approach for description of a propagation velocity has been 
proposed. It is strictly local and is based on the natural assumption of an ordinary 
measurement procedure. It does not need any a priori condition of periodicity, 
frequency, groups and packets or other canonical attributes. 

The definition elaborated is applicable in a natural way for arbitrary. In a 
mathematical sense it describes a propagation of perturbation in any geometrical 
field. Examples could be: an acoustic wave as a pressure perturbation, a grav- 
itation wave on a fluid surface, a spin wave in a solid state, etc. In particular 
the formalism is very suited for the description of particle propagation in field 
theory, where particles are considered as field perturbations. 

This approach results in the series of measurable propagation velocities. In 
zeroth order the propagation velocity coincides with the ordinary phase velocity 
and appears to be ordinary Lorentz covariant; further application gives rise to 
generate " Lorentz covariance of higher orders" . 

For propagation in inhomogeneous media, for instance for light in a medium 
with a space-dependent index of refraction phase velocities of higher orders ex- 
hibit an essential dependence on time component of the wave vector solely as a 
result of inhomogenity, treated canonically as a "dispersion". This appears to 
be a more general phenomena, and does not presuppose any periodic frequency 
and dispersive properties of media. It should be called for this reason "dynamic 
separation" . 
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